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We consider the Laplacian in a strip M. x (0, d) with the bound- 
ary condition which is Dirichlet except at the segment of a 
length 2a of one of the boundaries where it is switched to 
Neumann. This operator is known to have a non-empty and 
simple discrete spectrum for any a > 0. There is a sequence 
< a\ < a2 < ■ ■ ■ of critical values at which new eigenval- 
ues emerge from the continuum when the Neumann window 
expands. We find the asymptotic behavior of these eigenvalues 
around the thresholds showing that the gap is in the lead- 
ing order proportional to (a — a n ) 2 with an explicit coefficient 
expressed in terms of the corresponding threshold-energy res- 
onance eigenfunction. 
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1 Introduction 



Spectra of Dirichlet Laplacians in infinitely stretched regions such as 
planar strips or layers with local perturbations were studied recently in 
numerous papers. The motivation for this work came from applications 
in condensed matter physics, and also from the fact that it was itself 
an interesting mathematical problem. 

One of the simplest systems of this type is a free quantum particle 
confined to a pair of straight parallel strips with Dirichlet boundary 
conditions coupled laterally by a window in the common boundary. If 
they are of the same width d, one can employ the mirror symmetry and 
concentrate on the nontrivial part which is equivalent to the analysis 
of the Laplacian in a single strip with the Dirichlet boundary condition 
switched to Neumann at a finite segment of one of the boundaries. 

Such a system has at least one bound state for any "window" 
length 2a > as it was found in [ ESTV] and independently in ||BGRS . 
The discrete spectrum is simple, the eigenvalues A n , n = 0, 1, . . . , are 
continuously decreasing as functions of a and their number is linear 



in a up to an error term as the window is widening | ESTV|| . These 
properties follow from a simple bracketing argument which allows to 
squeeze the eigenvalues between those of a box covering the "coupled" 
part with Dirichlet and Neumann conditions at x\ = ±a. It shows, 
in particular, that there are critical values a n , n — 0, 1, . . . , at which 
new eigenvalues emerge from the continuum. By [ |ESTV| , |BGRS|1 we 
have ao = while generally we know only that 



/ nd 

V7i 



(n+l)d 



n 



1,2, 



This tells us nothing about the behavior of the eigenvalues around the 
critical points. 

The weak coupling asymptotics was studied for the ground state. 
A variational method of ||EV|| yields for small a a two-sided estimate 
of A (a) between two multiples of a 4 . This is indeed the leading term: 
using the matching method of [U, |Gal|| Popov derived in ||Po|| the 
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expansion 



While his argument is not fully rigorous because an estimate of the 
error term is missing, the formula itself raises no doubts, in partic- 
ular, because of its excellent agreement with the numerical result of 
ESTVfl. The result is subtle: recall that (|L2p differs substantially from 



the asymptotics corresponding to a local change in mixed boundary 
conditions, where the Birman-Schwinger technique is applicable and 
the leading term is a multiple of the window width squared ||EK|| . 

In the present paper we address the question about the behav- 
ior of the higher eigenvalues X n , n = 1,2,..., in the vicinity of the 
critical points a n . It is a natural counterpart of the coupling constant 
threshold problem of [[KS||. Our main result is that the gap between 



an eigenvalue \ n and the continuum is proportional to (a — a ri ) with 
a coefficient given explicitly in terms of the corresponding threshold- 
energy resonance eigenfunction. This fits well into the analogy between 
our problem and spectral properties of one- dimensional Schrodinger 
operators. The latter extends to higher dimensions, but the argument 
becomes more complicated and we leave that to another paper. 

To finish the introduction, let us say something about the method. 
We deal with a perturbative problem with respect to the parameter e 
defined as the excess of the window halfwidth a over the critical value 
a n . Since the latter is positive for n > 1, the problem in question is 
a regular one. This fact makes it possible to map the problem into 
an equivalent one with a small and local perturbation of the equation 
and the boundary condition fixed, i.e. independent of e. This is what 
we are going to do. We employ the technique introduced in ||Ga2 



for calculations of the eigenfunctions for one-dimensional perturbed 
Schrodinger operator; its advantage is that we arrive at the sought 
asymptotic formula in a straightforward and reasonably simple way. 
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2 The main result 



To formulate the result we need first to introduce some notation and 
recall a few simple facts about the problem in question. We employ 
Cartesian coordinates, x = (xi,X2), in which S = {x : < x 2 < d}. 
The upper strip boundary is T = { £ : x 2 — d}, while the lower 
decomposes into the union of r y a = { x : | x\ \ < a , x 2 = } and 
r a = { x : | xi\ > a , X2 = }. The operator H a we are going to 
consider is the Laplacian, H a ip = —Aip, in £ 2 (£) with the Dirichlet 
boundary condition on T a U T and Neumann on r y a \ it is a well-defined 
self-adjoint operator - cf. |Da| , Chap. 7], 



Proposition 2.1 The discrete part of a(H a ) is simple and the eigen- 
values A n (a) : (t^) < Ao(a) < Ai(a) < ■•• < (5) are continuous 
and monotonously decreasing with respect to a. There are numbers 
= a < ai < a 2 < ■ • • satisfying such that 

(a) for a G (a n -i) a n] the operator H a has exactly n eigenvalues, 

(b) for a > a n we denote e := a — a n , then the eigenfunction ipn^ 
associated with A n (a) has a definite parity with respect to x\, namely 

i/>V(-x 1 ,x 2 ) = (-l) n il>V{x 1 ,x 2 ), 

(c) for a = a n the equation —Ai[) = (f) 2 ^ with given boundary con- 
ditions has a solution G H^ 0C (T,), unique up to a multiplicative 
constant, which behaves like 

^\x) = Cl (±l) n sin (^) + O (e- 5 ^) (2.1) 

as x\ — > ±oo ; where 5 := ir\/3/d (in what follows we set c\ := a/2/7t 
for the sake of definiteness). 



PROOF: Most part was demonstrated in JESTV|, it remains to check 



the claim (c). Any solution can be expanded in terms of the transverse 
eigenfunction bases 

, x /2 . fnjx 2 \ [2 . f n(2j-l)(d-x 2 ) \ 
Xj (x 2 ) := ffln ^— j , Ux 2 ) := sm ^ j 
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with j = 1,2, ... . Since the problem has a mirror symmetry with 
respect to Xi = 0, it is sufficient to discuss the halfstrip part, X\ > 0, 
with the appropriate boundary condition at x\ = 0. Let is consider 
the even case. A solution of energy e (?) expresses as 

OO OO , / \ 

' cosh »,a 

3=1 3=1 J 

(2.2) 

for | Xi\ > a and \x\\ < a, respectively, where qj := \\/ j 2 — e andp,,- := 



f J(j - |) — e. The coefficients in the above relation are determined 
by the requirement of smoothness of ip at the segment X\ = a; we have 



bk (*i> ^ ' (*i> = o — t; 772 > ( 2 - 3 ) 

71 ^-(^-f) 



fc=i 



and 6 = {&,} is given as solution of an infinite system of equations 
which can be written concisely in the operator form 

Cb = (2.4) 

with 

C jk := [qj + p k tanh(p fc a)] (xj, <pk) ■ 

The odd case is similar, just cosh and tanh are replaced by sinh and 
coth, respectively. The allowed values of e are those for which a solution 
to the system ( |2.4| ) exists. 

We know from [ ESTV | that the sequence corresponding to an iso- 



lated eigenvalue of H a , i.e. e G Q, l) , belongs to £ 2 (j~ r ) for any r > 1, 
and that C = C(a, e) is Hilbert-Schmidt on C 2 {j~ r ) with r large enough 
independently of a and e. Choosing such an r it is straightforward to 
check that (a, e) i— > C(a, e) is jointly continuous in the corresponding 
Hilbert-Schmidt norm. Take a > a n and e n (a) := (-) A n (a), clearly 
e n(°) 1 as a -> a n +. The said continuity implies that the equa- 
tion (|2.4|) has for C(a n , 1) a unique solution in £ 2 (j~ r ), and by ( |2.3|) 
it determines a sequence b G ( 2 {j~ r )- Together they yield the sought 
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threshold resonance solution for a = a n with the asymptotics fl2.1| ) 
following from the first one of the relations (|2.2|). | 



Remark 2.2 The function tpn described in Proposition |2.1j has the 
following smoothness properties. It is infinitely differentiable every- 
where in £ except the endpoints of the segment 7 a ™. At these points 
the asymptotic formula 



(o), 



X 



l a n ri /2 sin^ + 0(r ± ) 



r± 







(2.5) 



is valid, where (r±, 9±) are polar coordinates associated with the vari- 
ables (±x\— a n ,X2) and a n is a some number (a unique one provided 
we fix ci). These asymptotics formulas can be verified in two easy 
steps. First, one should extend the function ifj^ to the mirrored strip 
{x : — d < X2 < 0} in the even way. This leads to solution of the 
equation — Aipn* 1 = (^) 2 ipn" 1 in the double strip {x : \x 2 \ < d} \ {x : 
\x\\ > a,x 2 = 0} with the Dirichlet boundary condition at the outer 
boundary and at the cut {x : \x%\ > a, x 2 = ±0}, and it is sufficient 



to employ the results established in [|NP|| for such elliptic problems 



Now we can state our main result: 

Theorem 2.3 The eigenvalue A n (a) of H a with n > 1 has the follow- 
ing asymptotic behavior, 



A n (a) 



H 2 n (a-a n ) 2 + O ((a-a n ) 3 ) 



as a — > a n +, where the coefficient is given by 



or alternatively by 



1 / 


M 0) 


a n J 


dxi 


s 






7ra 2 n 


/in : 


4 



dx\dx 2 , 



(2.6) 



(2.7) 



(2.8) 
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where a n is the number appearing in f\2.3[). The associated eigenf unc- 



tion tpn^ can be normalized in such a way that it satisfies the relation 

V>« = ^ + 0(e) , (2.9) 
in H 1 ((—R, R) x (0, d)) for any R > behaving asymptotically as 



^)( x ) = Cl (±l) n e- £ ^ |xi| sin ( 



71X 2 



0(e 



(2.10) 



when x\ — > ±oo, with S defined in Proposition pP 



Remark 2.4 The function ipn belongs, of course, to L 2 (S), but it 
does not have a limit in this space as e — > 0+ since the norm 11^4 ||l 2 (£) 
explodes in the limit. Furthermore, a n is nonzero for any n > 1. This 
fact can be easily deduced from the assertion (|2.7|) . Indeed, the as- 
sumption a n = implies immediately that ipn^ is independent on x\. 
However, this contradicts to the boundary value problems that ipn^ 
satisfies to. The coefficient a n being nonzero, the formula (|2.8| ) shows 
that the asymptotics of each A n is nontrivial in the leading order. 



3 Proof of Theorem 2.2 



The rest of the paper is devoted to the proof of Theorem |2.3| . With the 
scaling behavior of the spectrum in mind we can put d = n without 
loss of generality in the following. Furthermore, the mirror symmetry 
with respect to x\ — makes it possible to consider the halfstrip 
problem with the Dirichlet or Neumann condition at the cut. 

We need to introduce some notations. Let R\ = {x : x\ > 0} 
be the open right halfplane. As indicated above, we will work in the 
halfstrip II := EnR^., similarly we introduce 7 e := j an+£ n R\, T + : = 
T fl M+, and T £ := r a " +e fl R+. Moreover, we need a symbol for the 
cut at the symmetry axis, 7 := {x : X\ = , < x 2 < n}. Since we 
consider a fixed eigenvalue, we shall omit the index n when there is 
no danger of misunderstanding. 
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Let us outline the strategy of the proof. In the first step we are 
going to analyze the problem 



[A + l)u = m 2 u + f , are II, (3.1) 

l x u = , x e 7 , 
u — o , x g r u r + , 

du 



dx-. 



, x G 7o 



with a fixed function / at the right-hand side of the equation. The trace 
operator l x in the boundary conditions is defined as (l x u)(xi, x%) = 
«(0, X2) if n is odd and (l x u)(xi, x 2 ) = 5^7 (0, x 2 ) if n is even. The pa- 
rameter m is assumed to be complex and to lie in a (sufficiently small) 
neighborhood of zero (we indicate this neighborhood by T>). We also 
suppose that / is an arbitrary function from L 2 (J1) with a compact 
support. We will construct a solution of the problem (|3.1|) , meromor- 
phic with respect to m, with the following asymptotic behavior far 
from the cut of the half strip, 



dm) e~ mxi 



smx 2 + 0(e- V ^ Ixi ) , xi^+oo, (3.2) 



where c(m) is a constant determined by the function /. In the second 
step we will transform the original boundary value problem for the 
eigenfunction ip^ to another one with an equation the coefficients of 
which depend smoothly on e and the boundary condition is indepen- 
dent of e. What is important is that the reformulated problem will be 
of the form fl3.1|) with a particular right-hand side f — f e for which we 
deduce a sufficient explicit representation. Combining the latter with 
properties of the solution to ( |3.1| ) will finally obtain the announced 
results concerning A e and tp( £ \ Since the proof of Theorem |2.3| divides 
in this way naturally into two steps, we shall discuss them separately 
in the two following subsections. 



3.1 Solution of the problem (|3.1|) 

As we have indicated our aim is to construct a solution of (O) which 
is meromorphic in m G T>. Let us say more explicitly what we mean 
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by that. It is easy to see that there is a unique solution for m G 
V R {m : Rem > 0} which decays as X\ — > oo. We shall check that as 
a function of m it is analytic in T> n {m : Rem > 0} and extend it to 
the remaining part of T>. The extension will be for Rem < again a 
solution to (|3.1|) with the asymptotics (|3^ ), and in addition, it will be 
meromorphic mm ET> with just one simple pole at m = 0. Of course, 
the extension will be bounded at large distances if Re m = and will 
be increasing for Rem < 0. Speaking about solutions everywhere in 
the following we mean always such analytic continuations. 

Recall that a function F with values in some Banach space X is 
said to be holomorphic in T> if it is different iable (in the sense of the 
norm of X) at each point m G T>. If this function is holomorphic in 
T> everywhere except a discrete set of points which are poles of F, i.e. 
Laurent's series of F at such a point has at most a finite number of 
negative terms, then F is said to be meromorphic. 

We introduce some notations. Let C(X, Y) be the Banach space 
of bounded linear operators from a Banach space X into a Banach 
space Y. We will use the symbol TC(X) for the class of functions with 
values in X which are holomorphic with respect to m G T>, and M.(X) 
for the class of meromorphic on m G T> functions with values in X. 
For the sake of brevity we also introduce the notations 7i(X,Y) := 



||SP| , Sec. XVI. 4]. Let R be a fixed number larger than a n , Ur := 
LT n {x : x\ < R}, and g a function from L 2 (I1r) which can be also 
regarded as an element of £ 2 (II) if we put it equal to zero for x\ > R. 
The problem 

(A + l)v = m 2 v + g , xeU, v = 0, xG<9IT, (3.3) 
can be easily solved by separation of variables, the solution being 




m ) = - HTT / G k (x,t,m)g(t)d 2 t, 
k=i k i 




(3.4) 



G k (x,t,m) = (e~ Mklxi - tl1 - e - Mk{xi+tl) ) sin fct 2 sin fcx; 



2 • 
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where Mi = m and M k = \/k 2 — 1 + m 2 for > 2. For the sake of 
brevity we use the notation d 2 t = dt\ dt 2 . Obviously, the formula ( |3.4|) 
is valid for all m G T>, not only for Rem > 0. The function v can 
be represented as v = T 1 (m)g, where Ti{m) : L 2 (Jl R ) — > H 2 (J1^) is 
a bounded linear operator for any positive R. It is straightforward to 
check that v G n(H 2 {Ii^)) and r a (-) G H(L 2 (n fl ), # 2 (n 5 )) in V for 
any positive R. 

In the next step we consider another boundary value problem for 
an unknown function w, namely 

Aw = Av, xeU R , (3.5) 

w = o, x g (r + u To) n du R , 

w = v , Xi = R, 
dw 

— = , x G 70 , 

G>£ 2 

/ x w = , x G 7 . 

Since t> G -ff 2 (II^) we have At> G L 2 (II/ ? ), and thus the problem 
(3.5) has a unique solution w G ^(n^) - see, e.g., fLa] , Sec. 2.5, 



Rem. 5.1]. Using the standard theorem on smoothness of solutions of 
elliptic boundary value problems we can conclude that w G H 2 {Tl RtS ) 
holds for each s > 0, where Il^ s := U R \D S with D s := {x : x 2 > , 
(zci — a n ) 2 + ^1 < s<2 }- Hence we can define the linear operator T 2 by 
w =: T 2 v which is a linear bounded map from H 2 (U R ) into H 1 (H R ), 
and from H 2 (U R ) into i7 2 (il^ s ) for each fixed s. 

Let Xr{ x i) be an infinitely differentiable mollifier function equal 
to one for x\ < R — 1 and vanishing for xi > R. We take the two 
functions considered above and construct u as a smooth interpolation 
between them, 

u(x, m) := xr{x\)w(x, m) + (1-xr(xi))v(x, m) . (3.6) 

Since w(x,m) = T 2 T 1 {m)g, we have w G HiH^ILn)) n H(i/ 2 (n fliS )) 
as a function of m for each s > 0. Thus we can introduce the operator 
Ts(rn) which maps a function g G L 2 (U R ) to the function u deter- 
mined by (|3.6D, where v and w are the solutions of (^.3|) and 
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respectively, for this g. It easy to see that T 3 (-) G H(L 2 (Il R ) , H l {U.^)) 
and T 3 (-) e 7^(L 2 (IIr), if 2 (n^ g )) in D for any pair of positive R, s. 

According to the definition of v and w, the function u satisfies all 
the boundary condition of ( |3.1| ). Applying the operator (A + 1 — m 2 ) 
to this function, we obtain 

(A + 1 - m 2 )u = g + (w-v)(A + 1 - m 2 )xn + 
+ 2 (V x xr, ^x(w-v)) r2 , 

where we have used in the calculation the equations which v and w 
satisfy This result shows that the function u solves the problem ( |3.1| ) 
if and only if g satisfies the following equation, 

g + T 4 (m)g = f, (3.7) 

where 

T 4 (m)g := (w-v)(A + 1 - m 2 )xR + 2 (V^Xr, V x (w-d)) r2 ; 

recall that both w and v are obtained from g by actions of the operators 
specified above. 

Proposition 3.1 The operator T^m) is compact for any R > as an 
element of £(L 2 (U R ), L 2 (ITr)) and the function m \— > T 4 (m) belongs 
to H(L 2 (Tl R ), L 2 (Tl R )) inV. 

Proof: We denote 

T 41 (m)g :=(w-v)(A + l-m 2 )x R , 

T 42 (m)g :=2 (V^h, V x (w-f )) Ra + (w-w)A^ . 

Using the described properties of Tx(m) and T 2 it is easy to see 
that Tn{m) is a bounded linear map from L 2 (U R ) into /^(II/j). The 
operator-valued function T 4 i(-) belongs to 7i(L 2 (ITH), H l (JX R )). The 
function is smooth and its support lies in IlR\nR_i. Using this to- 
gether with the fact that the function w—v belongs to if 2 (n J j\ri/j_i), 
we conclude that T 42 (m) : L 2 (ITr) — > ^(Hr) is a linear bounded 
operator belonging to 7^(L 2 (ITr), ^(IIr)) as a function of m. Hence 
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r 4 (.) = T 41 (-) +T43O) G 7^(L 2 (IIr), ^(IIr)), and therefore T 4 (m) 
is compact when considered as an operator from L 2 (ITr) to L 2 (Il/j), 
belonging to W(L 2 (n fl ), £ 2 (n fi )) w.r.t. m. I 



Proposition |0| shows that the equation ( |3.7| ) can be studied using 



Fredholm theorems. This will help us to solve our original problem: 
to construct the solution of ( |3.1| ) we have to solve the equation (|3.7| ), 



then by the procedure described above its solution gives rise to the 
solution u of the boundary value problem ( |3.1| ): u = T 3 (m)g. 



Proposition 3.2 The problems (3.7) and (3.1) are equivalent: to each 



solution g of \3. % ) there is a unique solution u = T 3 (m)g of ( \3.1\ ), and 
vice versa, for each solution of (\3.1\) there exists a unique g solving 



( \3. Tj) such that u = T 3 (m)g. The equivalence holds for any m G V. 



Proof: The first part has been proved above, it remains to check 
invertibility of the operator T3. Let u be a solution of the problem 
(3~I). Notice that the function u is infinitely differentiable outside the 



support of / since there it is a solution of a homogeneous equation. 
We have to construct the solution g of (|3.7|) such that u = T 3 g. Let 
us first determine the functions v and w. By U we denote a solution 
of the following problem, 

AU = 0, xeU R , (3.8) 
U = 0, xedU R \^U^, 
U = u, XG70U7, 

which is unique in if 1 (n fi ), and moreover, it belongs to C°°(Il R ) by 
a standard result on the smoothness of solutions of elliptic equations. 
We set 

v(x, m) := u(x, m) — xr(%i)U(x, m) , 
w(x, m) := u(x,m) + (1 — Xr{xi))U(x, m) . 

One checks easily that the functions u = xr w + (1 — Xr) v i w an d v 
satisfy all the required boundary condition and that Aw = Av holds 
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in U R . Now it suffices to use the equation of the problem ( |3.3| ) to 
determine the function g by 

g(x, m) :— (A + 1 — m 2 )v(x, m) 

= (A + 1 — m 2 )(u(x, m) — xr(xi)U(x, m)) 
= f(x) - (A + 1 - m 2 ){ X R(xi)U(x, m)) . 

Let us check that this g solves ( |3.7| ). Using the definition of U, v, and 
w, we compute directly the action of T^m) on g obtaining 

T A {m)g — U(A + 1 — m 2 ) XR + 2{V xX r, V x U) 

= U(A + 1 - m 2 ) XR + 2(V xX r, V X U) + Xr AU 
= (A + l - m 2 )( XR U). 

The last two relations show that g is a solution of the equation (|3.7|) . 
Let us check the uniqueness: suppose that there are two solutions 
gi and g% of (|3.7| ) leading to the same function = u^ 2 \ Then 
g := g\—g% ^ gives rise to u := u^—u^ = 0. Let v and w be the 
solutions of (|3~3| ) and ( |3.5| ), respectively, associated with u, and put 
U := w — v. Then it easy to see that this U solves ( fj.8[ ), where the 
boundary condition on 7 U 7 is homogeneous. Such a solution of ( |3.8D 
is unique, U — 0. Thus to = t> = it = holds in Ur, and therefore 
<7 = (A + 1 — m 2 )v = in which is a contradiction. | 

The solution of eq. ( |3.7|) depends on m, i.e. g — g(x,m). Our next 
aim is to clarify a nature of this dependence and to look what this 



implies for the solution of the problem (3T). We employ the following 



result borrowed from ||SP| , Sec. XV. 7]: 



Theorem 3.3 Let T> be an open connected domain of the complex 
plane of the variable m and {T(m) : m G T>} be a bounded holo- 
morphic family of compact operators from the Banach space X into 
itself. Moreover, assume that there exists a point m* G D such that 
(I + T{m^))^ 1 G C(X,X). Then m 1— > (I + T(m))~ 1 is a meromorphic 
function in T> with values in C(X,X). 



In our case the Banach space mentioned in Theorem |3.3| is L 2 (I1^) 
and T{m) = T 4 (m). The existence of m* is easy to establish: since the 



13 



equation ( |3.7| ) is equivalent to the boundary-value problem by 
Proposition |3.2| , it is sufficient to prove the existence of m* for which 
the problem ( |3.1[ ) has no nontrivial solution for / = and m = m* 
with the asymptotics ( |3.2| ). This is true for m* > which is sufficiently 
small, because assuming the contrary would lead us to the conclusion 
that H an has the eigenvalue (1 — m%), and this in turn would contradict 



to the claims (a) and (c) of Proposition |27i| . The compactness of Tjjn) 
and its holomorphic dependence on m follow from Proposition ft.lj 

Thus we may apply Theorem [3.3| to the equation (|3.7p . We implies 
that (J + T^{m))~ 1 exists and it is meromorphic as an operator- valued 
function in V, (I + T 4 (-))" 1 G M(L 2 (Tl R ), L 2 (Tl R )). Poles of this 
function are the values of m for which the equation ( |3.7| ) with the 
vanishing right-hand side has a nontrivial solution. The value mo := 
has this property as it follows from the claim (c) of Proposition |2.1| . 
Let 0o be a solution of the equation (|3.7|) for / = and m = mo. The 
function O is unique up to a multiplicative constant, the uniqueness 
being implied by Proposition |2J](c) and Proposition p.lj The remain- 
ing ambiguity is removed if we set ^ < -°- ) = T 3 (O)0o- Making the domain 
T> smaller if necessary we can achieve that zero is the only pole of 
the function (/ + T^-)) -1 contained in T>. In such a case the solution 
g = (I + T4(m)) _1 of the equation fl3.7p can be for any nonzero m ET> 
represented as 

9 = g{m) = ^=f +g(m) , (3.9) 
vnr 

where the integer k > 1 is the order of the pole, the functions g(m) 
belong to L 2 (ITr) for nonzero m G T>, and g(-) G A4 (L 2 (I1r)) may have 
a pole at zero of order not exceeding k — 1. 



Let us stress that Theorem says nothing about the order of 



this pole. Next we are going to prove that the pole in ( |3.9| ) is simple, 
i.e. k = 1. Substituting (|3.9| ) into ( |3.7|) and comparing the leading 
terms in m we see that = K [f](f) , where K [f] is a constant 
depending on /. The representation ( |3.9|) in turn yields the following 
expression for u: 

u = T 3 (m)g = r hum . (3.10) 

m K 
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Here u(-) is a function belonging to A^if^IIjj)) fl M (H 2 (Il^ J) for 

any positive i?, s which again may have pole at zero of order not 
exceeding k—1. Multiplying the equation in the problem ( |3.1| ) by ip^ 
and integrating by parts we see that 



f^d 2 x = -m 2 I ^ud 2 x + I ( V (0) ^ - 



u- 



dxi dx 



II £ 



dx 2 . 

xi=R 

(3.11) 



For <i sufficiently large we have 

oo oo 

U \ X1 =R = J2 C j( m ) ^jx 2 , ^ (0) L 1= r = 53 C J sin jx 2 

oo 

53 \/j 2 — 1 +m 2 Cj (m) sin jx 2 , 

oo 



9m 



0^,(0) 



dx i 



xi=_R 



j=2 



where the functions Cj in view of ( p.lOj ) behave as 



Combining the above relations and using the normalization condition 
Ci = \plpK we deduce that 



{0) du 

^ a u ^ 

ax\ ox i I 



xi=R 



dx, = -^Q + 0(m- k+2 ) 



The first integral at the right hand of ( |3.11| ) behaves as 



m 2 / ^ {0) ud 2 x 



m 



k-2 



\^\ 2 d 2 x + 0(m- k+3 ). 



II £ 
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in the limit m — > 0. Substituting the expressions obtained above into 
(|3.11| ) and comparing the coefficients at the powers of m, we conclude 



that k — 1, and furthermore, that 

K [f] = - [ M (0) d 2 x = - I d 2 x, (3.12) 



To get the last relation we have used the fact that / has a compact 
support by assumption which therefore lies in IT^ for R sufficient large. 

Next we denote A(m) = T 3 (m)(I + T 4 (m)) _1 . It follows from the 
relations ( |3TT0|) , ( |3TT2D that 



u = A(m)f =: A (m)f + A 1 (m)f = K [f}^- + A 1 (m)f, 

Ai(-) e n(L 2 (u R ),n(H\u n )) n n(L 2 (u R ), h 2 (u Rs )), 

where R and s are arbitrary positive. 

It is convenient to summarize all the conclusions made above in 
a single theorem which represents the main result of this subsection. 

Theorem 3.4 Let f 6 L 2 (U R ), and assume that R, s are arbitrary 
positive numbers, m GD. Then for m G T> \ {0} there exists a unique 
solution of the boundary value problem (\3.^ ) given by u = A(m)f. As 
a function on m this solution belongs to 7W(if 1 (n^)) H M.(H 2 (U^ s )). 
The neighborhood T> can be chosen in such a way that the function 
A(-) G M(L 2 (n R ),H\n n ))nM{L 2 (U R ),H 2 {Ii n ;)) has just one pole 
m = which is simple. The operator A{m) can be decomposed into the 
sum A(m) = Ao(m) + Ai(m), where the operators Ai(m) are defined 
by ( \3.1c\ ). Finally, the relation ( \3.13j ) holds true. 



3.2 The asymptotic analysis 

To find the behavior of the quantity A n (a) around the threshold we 
will analyze m £ defined as m 2 := 1— X n (a n + e), which satisfies m £ — > 
as e — ► 0. Let y = (3/1,3/2) be the Cartesian coordinates of a point in 
the strip. The boundary value problem problem for the eigenfunction 
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can be then written as 

{A y + 1)^ (£) = m 2 e i) {e) , yen, (3.14) 

^ = o, yer £ ur + , 

— — = , 1/6%. 
dy 2 

As we have announced, we want to get rid of the ^-dependent bound- 
ary condition passing to one independent on e, while the correspond- 
ing equation will have coefficients which depend smoothly on e. To 
construct the appropriate transformation we use an infinitely differen- 
tiable mollifier function \ of the variable y\ which is equal to one for 
Vi ^ IP2, P3) an d vanishes for y x < Pi an d Vi > Pa- Here Pi, i — 1, . . . , 4, 
are positive constants, Pi < P2 < P3 < Pa, such that P2 < a n +e < P3 
holds for all e from a fixed neighborhood of zero. We consider the 
function x\ : R + — > R + defined as 

x x : x^y^e) = y x -ex(yi) ■ 

It is easy to see, in particular, that x± = a n as y\ = a n + e. Taking the 
first two derivatives we get 



dx\ t , . d X\ n 

dyi dyf 



from where it follows that for all sufficiently small e the first derivative 
of x\ w.r.t. yi is nonzero everywhere in 1R + . Thus the map x\ is bijective 
and can be used to define a change of the variable, y\ 1— > £1(2/1, e). 
The problem fl3.14| ) expressed in the variables x = (xi, X2), X2 '■— y2, 
becomes 

(A x + eL £ + l)^ £) = m 2 £ ^ £ \ xen, (3.15) 

l x ^ e) = , x e 7 , 

= , x e r u r+ , 



dxo 



, x e 70 . 
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The operator L £ appearing in the transformed equation is defined by 

.— -2 +6i(a?i,e) — 

OX\ OX\ 

bxx{x x ,e) := -2y^(y 1 (x u e)) + e(x / (yi(»i, ^))) 2 
6i(o;i,e) := e)) • 



6n(a;i,e)- 5 -2 +6i(a?i,e)- 5:r , (3.16) 



The functions fen and h\ are obviously infinitely differentiable, they 
vanish for y\ ^ [fli, P2) U [/3 3 , /5 4 ] and satisfy in the limit e — > the 
asymptotic formulae 

buixue) = -2^{ Xl ) + 0(e), 6i(xi,e) = + O(e) 

uniformly in the variable 2^. 

Now we can proceed to the calculation of m £ and ip( £ \ Propo- 



sition |2.1| ensures that the eigenfunction and eigenvalue exist. The 
function ip^ decays as \xi\ — > 00 and A n (a n + e) is real. These two 
facts imply that m £ is real and positive. It is obvious that ip^ solves 
the problem for m = m £ and / = —£L E ip^ £ '. This allows us to 
seek for the solution of ( |3.15|) in the form 

^ = A(m £ )f £ , (3.17) 

where the function f e is assumed to an unknown element of L 2 (ITr) 
with R > (3&. This may appear strange at a glance, because we know 
that f £ = —eL £ ip( e \ however, we want to obtain another formula for f £ 
not involving ip( £ \ The function ip^ defined by ( p. 17 ) satisfies all the 



boundary condition of the problem ( |3.15| ). In order to be a solution 



of (|3.15|) , it is necessary and sufficient for tp^ to be a solution of the 
corresponding equation. Substituting into the latter the formula for 
we arrive at the equation for f £ which reads 

(I + eL £ A(m £ ))f £ = 0. (3.18) 

In view of (|3.13| ) we have 

A(m £ )f £ = K [f £ ]^— + A l {m £ )f £ . (3.19) 
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We substitute this representation for A{m e )f £ into ( |3.18| ) obtaining 
f e + —K [f £ ]L £ ^ + eL £ A x {m £ )f £ = . (3.20) 



It is clear that L e A\(-) G H(L 2 (I1 R ), L 2 (I1 R )) and as a function of 
(m, e) the operator L e Ai(m) is jointly continuous. Thus for sufficiently 
small e the inverse operator B(m,e) = (I + eL £ Ai{m))~ l : L 2 (Jl R ) — > 
L 2 (Ur) exists and converges to the identity map as e — > uniformly 
in m. It is also obvious that B(-,e) G H,(L 2 (J1 R ), L 2 (I1 R )). Applying 
now B(m £ ,e) to the equation (|3.20|) we find 

f e + —K [f £ ]B(m £ , e)L £ ^ = . (3.21) 

Tfl £ 



Applying further K to ( |3.21| ), we get one more equation, 

K [f £ ] + —K [f £ }K [B(m £ ,e)L £ ^} = 0. 

m £ 

Notice that K Q [f £ ] can not be zero, because otherwise ( p.21| ) would 
imply f £ = which yields ip^ = 0. The last equation induces the 
relation 

m £ = -eK [B(m £ , e)L £ ^} , (3.22) 

which can be regarded as an equation for m £ . It is easy to see that 
B(-,e)L e ipW G H{L 2 {U R )) and the function (m,e) i-> B(m 1 e)L e ipW 
is jointly continuous. This immediately implies that K [B(m, s)L £ ip^] 
is a holomorphic with respect to m and jointly continuous as a function 
of (m, e) . Consequently, for all e small enough the estimate 

e \[K [B(m,e)L £ ip {0) }\ < \m\ , m G dV , 

holds true. Using this inequality in combination with the Rouche theo- 
rem we conclude that the functions h\ : hi(m) = m and h,2 : h2(m) = 
m + [Ko[B(m,e)L £ tp^] have the same numbers of zeroes inside T>. 
This means that the equation (|3.22| ) has a unique solution in T> for all 
sufficient small e. On the other hand, due to Proposition |2.1| we know 
that there is a root of equation ( |3.22| ) that converges to zero, namely 

m £ = ^J\ n (a n +e) - 1 > . 
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Consequently, it is the only root of ( |3.22| ). Thus the function 

f e = -eB{m e ,e)L e 1>®, (3.23) 

where m e is the solution of ( |3.22j ), solves the equation ( |3.18| ). It means 
that the function f e defined by ( |3.23|) gives rise to the eigenfunction 
tp( e ) = A(m e )f £ corresponding to the eigenvalue X n (a n + e) = 1 — m 2 . 
In fact, we have also proved that there is just one value of m = m £ 
tending to zero as e — > 0, for which the boundary value problem ( |3.14|) 
has a nontrivial solution. This solution decays as \x\\ — > oo, i.e. there 
are no non-decaying or even increasing solutions. 

The equation ( p.22[ ) allows us to calculate the asymptotic expan- 
sion for m £ . Since B(m, e)L e ip^ is holomorphic with respect to m and 
jointly continuous in (m, e), and since m £ tends to zero as e — > 0, then 
by the equation ( |3.22| ) we have 



m e = en + 0(em e ) , e — > 
L = -2 X '(x 1 )-^-x"(x 1 ] 







Combining these relations, we can rewrite the error term as follows, 



m e = e/j, + 0(e 2 ) , e — 
Next we want to express the coefficient fi: 

d 2 . ,„ , d 



0. 



(0) 



2 *'^dxi+x"^ dxi 



^ d 2 



X . 



(3.24) 



(3.25) 



Using the equation which the function satisfies we find 



d 2 // 9 

2 x{xi)^j2 +x"(yi)- 



(0) 



dx\ dx\ 

Integrating then by parts and taking into account properties of ip^ 
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together with the definition of the mollifier x we have 



fi = -\im I ^ (0) (A + 1) I x( x i) 



n\D s 



d 2 x 



lim 

s 



J \ drdxi dx\ dr J 

dD s 



ds . 



In order to evaluate the last integral along dD s we replace ip^ by its 
asymptotics (|2.5|) and pass to limit s — * obtaining 



< r / sin 2 — d9 



which proves (|2.8|) with /i = \i n . In the same way we get 



II 



1 dx\ 



■(A + l)Vi 0) rf 2 x = a n ^- + 2 / ^^-^d 2 * 



/ 


M 0) 


n 





g? 2 x 



which yields the representation ( |2.7| ) for //, 



2 

(Xn. 



(o) 



9xi 



d x = — 



8^n 



(0) 



9xi 



In the last relation we employed the fact that ipn^ has a definite parity. 

Let us finally pass to discussion of the eigenf unction. In view of 
J319D , ( B ) we find that is equal to 



A{m £ )f £ 



(o) 
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Using now the equation ( |3.22j ) together with the fact that the function 
A 1 (m)B(m,e)L £ ip^ is holomorphic and continuous as a function of 
the respective variables, we conclude that the relation 



(o) 



^) = ^(°) _ e A 1 (m e )B(m e , e)L e ip 

is valid in H 1 (Un) for each R. This yields the relation fl2.9|) and the 
asymptotic behavior fl2.10| ) concluding thus the proof of Theorem |2.3| . 
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